Bright solitary waves in a Bose-Einstein condensate and their interactions 
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We examine the dynamics of two bright solitary waves with a negative nonlinear term. The 
observed repulsion between two solitary waves - when these are in an antisymmetric combination 
- is attributed to conservation laws. Slight breaking of parity, in combination with weak relaxation 
of energy, leads the two solitary waves to merge. The effective repulsion between solitary waves 
requires certain nearly ideal conditions and is thus fragile. 

PACS numbers: 03.75.-b, 03.75.Lm 



I. INTRODUCTION 

One of the many interesting features of Bose-Einstein 
condensed atoms is that they can support solitary waves, 
in particular when these are confined in elongated traps. 
Under typical conditions, these gases are very dilute and 
are described by the familiar Gross-Pitaevskii equation, a 
nonlinear Schrodinger equation with an additional term 
to describe the external trapping potential. 

It is well known that the nonlinear Schrodinger equa- 
tion (with no external potential) supports solitonic solu- 
tions through the interplay between the nonlinear term 
and dispersion. In the presence of an external trapping 
potential, the Gross-Pitaevskii equation becomes non- 
integrable. In elongated quasi-one-dimensional traps, it 
is reasonable to approximate the three-dimensional solu- 
tion of the Gross-Pitaevskii equation by separating lon- 
gitudinal and transverse degrees of freedom pj. The re- 
sulting effective one-dimensional nonlinear equation has 
a nonlinear term that is not necessarily quadratic [l], 0] • 
Still, such nonlinear equations support solitary-wave so- 
lutions, which must be found numerically. 

Solitary waves have been created and observed in 
trapped gases of atoms [1, 0, H, @ . In the initial exper- 
iments [jj IH the effective interaction between the atoms 
was repulsive. In this case, the solitary waves are lo- 
calized depressions in the density, which are known as 
"grey" solitary waves. These waves move with a velocity 
less than the speed of sound. When the minimum of the 
density (at the center of the wave) becomes zero, they do 
not move at all and thus become "dark" . 

More recently the two experiments of Refs. 0, 0| con- 
sidered the case of an effective attraction between the 
atoms and observed "bright" solitary waves, i.e., blob(s) 
of atoms which preserve their shape and distinct iden- 
tity. Strecker et al. @ created an initial state of many 
separate solitary waves. While these independent waves 
were seen to oscillate in the weak harmonic potential in 
the longitudinal direction, they did not merge to form 
one solitary wave. In other words, they behaved as if 
the effective interaction between two of these waves were 
repulsive. 

Numerous theoretical studies have been motivated by 



the experiments of Refs. [Ell], see e.g., Refs. @,H,@- Ref- 
erence 0] offered an explanation for the observed effective 
repulsion between solitary waves. As argued there, the 
experiments had been performed in a manner that gave 
rise to a phase difference of tt between adjacent solitary 
waves. According to an older study [l(|, solitary waves 
with a phase difference equal to tt indeed repel each other. 



In the present study we use a toroidal trap 11] as a 
model for examining the time evolution of a system that 
initially has two solitary waves using numerical solutions 
to the corresponding time-dependent one-dimensional 
Gross-Pitaevskii equation. Remarkably, such toroidal 
traps have been designed (T^ |. and very recently persis- 
tent currents have been created and observed in such 
traps [l3[ . The basic conclusion of our study is that the 
effective repulsion between solitary waves is due to con- 
servation laws and thus fragile. 



In what follows, we first present our model in Sec II. 
In Sec. Ill we examine the dynamics of the gas in the 
case of weak dissipation, starting with perfectly symmet- 
ric/antisymmetric initial conditions and with no external 
potential along the torus. We observe that the symmetric 
configuration of two blobs merges on a short time scale; 
the blobs in the initially antisymmetric configuration re- 
main distinct and separated. Using these results as "ref- 
erence" plots, we examine in Sec. IV the effect of a weak 
random potential on perfectly symmetric/antisymmetric 
initial conditions. We also examine in Sec. V the time 
evolution of states that deviate slightly from perfect sym- 
metry/antisymmetry in the absence of any random po- 
tential. In both cases, the symmetric (or nearly symmet- 
ric) initial configuration shows essentially the same be- 
havior as the reference symmetric system. On the other 
hand, the antisymmetric configuration with the addition 
of an extra weak random potential and the nearly anti- 
symmetric configuration with no external potential both 
lead to a merger of the two blobs after a moderate tran- 
sient time. In the antisymmetric case, the final state is 
strongly influenced by weak deviations from the "ideal" 
case. Finally, in Sec. VI we discuss our results. 
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FIG. 1: Plot of |#(0,i = 0)| and of <f>(6, t = 0) for the symmet- 
ric (higher) and for the antisymmetric (lower) case, a = ±1 
in Eq. ©. 
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FIG. 2: Snapshots of \9(6, t)\ (higher curve in each panel) and 
of 4>(6, t) (lower curve in each panel) of the order parameter 
* = \^!{6,t)\e l4 ' { - e ' t \ for the symmetric initial configuration, 
a = 1, for t/to = 0, 10, 50, 100, and 150. The axes are the 
same as in Fig. 1. In all the above graphs there is no external 
potential, V = 0. 
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II. MODEL 

We consider a tight toroidal trap and use the mean- 
field approximation. Tight confinement along the cross 
section of the torus allows us to assume that the trans- 
verse degrees of freedom are frozen, and thus the corre- 
sponding time-dependent order parameter 9(9,t) satis- 
fies the (one-dimensional) equation 

89 h 2 d 2 9 

where g — 8nNaR/ S , and V(9) is the external poten- 
tial measured in units of Eq = ft 2 / (2 MR 2 ). Here, M is 
the atomic mass, R is the radius of the torus, N is the 
atom number, a is the scattering length (which is taken 
to be negative), and S is the cross section of the torus. 
The total length of the torus is chosen to be 167T in our 
simulations. 

As shown in Refs. [3, El, below a critical (negative) 
value of the parameter g, there is an instability from 
a state of homogeneous density to a state with local- 
ized density that breaks the rotational invariance of the 
Hamiltonian. This localized state corresponds to a soli- 
tary wave, and the critical value of g is g c = —n for 
the parameters chosen here. We adopt the value of 



FIG. 3: The energy of the system in units of Eo = 
h 2 /{2MR 2 ), versus time for < t/to < 150 corresponding 
to Fig.H 

g = —4 < g c in all simulations below. Therefore, the 
lowest energy state is a single solitary wave. 

We add an extra term on the left side of the above 
equation to model dissipation and write 

The real positive dimensionless parameter 7 describes the 
"strength" of dissipation. Since we solve an initial value 
problem, we also need to specify the initial condition. 
This is 

9(9, t = Q)= J—x W - 9 ) + a^(9 + 9 )), (3) 
Vl + a 

Here ip(0) = A/cosh(A#), with A = 3/2, is a static, well 
localized blob. We choose 9q = 2tt/5 so that the two 
blobs are reasonably distinct but still have a small over- 
lap, as shown in the graphs of Fig.[T] for a = ±1. 

Figures | | 1 | H and H show snapshots of 
\9(9,t = 0)| and of <j>(9, t = 0) for the order parame- 
ter 9(9, t) = \^(9,t)\e^ e ^>. Figs.El El H HTJ and[n] 
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FIG. 4: Snapshots of 1^(0, t) | (higher curve in each panel) and 
of <f>(9, t) (lower curve in each panel) of the order parameter 
* = |\f r (0,f)|e i '* (9, * ) , for the symmetric initial condition, for 
t/to ~ 0,10,50,100, and 150, for a weak random potential 
(shown in Fig. ll4[) , with a symmetric initial configuration, a = 
1. The axes are the same as in Fig. 1. 
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FIG. 5: The energy of the system in units of Eq — 
h 2 /(2MR 2 ), versus time for < t/to < 150 corresponding 
to Fig.H 

show the corresponding energy of the gas as function of 
time. 



III. TIME EVOLUTION OF THE "IDEAL" 
SITUATION 

To understand the effects of a weak random potential 
and the effect of slight asymmetries in the initial condi- 
tion (to be considered in Sees. IV and V), it is instructive 
to start with the situation where there is no external po- 
tential, V{9) = and an initial configuration which is 
either perfectly symmetric (a — 1), or perfectly antisym- 
metric (a — —1), i.e., 

*(0, t = Q) = -^=[i>{9 - ) ± 1>{9 + 9 Q )]. (4) 

We also fix the value of the dissipative parameter equal 
to 7 = 0.05. Figures O and [5] show snapshots of 
\^(9, t)\ as well as the phase 4>(9,t) of the order pa- 
rameter ^(9, t) for the symmetric and the antisymmetric 
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FIG. 6: Snapshots of \^(0, t)\ (higher curve in each panel) and 
of 4>(6, i) (lower curve in each panel) of the order parameter 
* = |*(t9,t)|e i0(9 ' t) , for the almost symmetric initial condi- 
tion, a = 1.01, for t/t = 0,10,50,100, and 150. The axes 
are the same as in Fig. 1. In all the above graphs there is no 
external potential, V = 0. 
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FIG. 7: The energy of the system in units of Eo = 
versus time for < t/to < 150 corresponding 

toFig.H 



case, respectively. The snapshots shown in Fig. [2] cor- 
respond to t/to = 0, 10, 50, 100, and 150, and in FigM 
to t/to = 0,10,100,300, and 400. Here t = E Q /h = 
[h/(2MR 2 )}- 1 is the unit of time. Figures [3] and [9] 
show the energy of the system as function of time for 
< t/to < 150, and for < t/t < 400, respectively. 

As seen in these graphs, the symmetric configuration 
(Fig. ^ merges quickly into one soliton and as time in- 
creases, it eventually approaches the equilibrium solu- 
tion. On the other hand, the two blobs do not merge 
in the antisymmetric case (Fig.[SJ). This is a direct con- 
sequence of the fact that the initial configuration has a 
node at 9 = 0. Because of the symmetry between 9 and 
—9, = 0, t) must be zero for all times t > 0. As a re- 
sult, the two blobs never merge as a simple consequence 
of parity conservation. The parity operator commutes 
with the Hamiltonian, and parity is therefore a conserved 
quantity. Only numerical errors can eventually lead to a 
single-soliton profile (with lower energy). That this does 
not happen provides a check on the accuracy of our nu- 
merics. 
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FIG. 8: Snapshots of 1^(0, t) | (higher curve in each panel) and 
of <f>(9, t) (lower curve in each panel) of the order parameter 
* = |*(6», t)\e l4,{ - e ' l \ for the antisymmetric initial configura- 
tion, a = -1, for t/t = 0,10,100,300, and 400. The axes 
are the same as in Fig. 1. In all the above graphs there is no 
external potential, V — 0. 
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FIG. 9: The energy of the system in units of Eq = 
h 2 /(2MR 2 ), versus time for < t/to < 400 corresponding 
to Fig.|U 
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FIG. 10: Snapshots of 1^(6, t)\ (higher curve in each panel) 
and of <j}(8, t) (lower curve in each panel) of the order param- 
eter \t = |\E r (0,i)|e 1 * (fl, * ) , for the antisymmetric initial condi- 
tion, for t/to ~ 0, 10, 100, 300, and 400, for a weak random 
potential (shown in Fig.ll4[l. with an antisymmetric initial 
condition, a = —1. The axes are the same as in Fig. 1. 
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FIG. 11: The energy of the system in units of Eq = 
h 2 /(2MR 2 ), versus time for < t/to < 400 corresponding 
to Fig.[l0l 



IV. EFFECT OF THE RANDOM POTENTIAL 
ON THE TIME EVOLUTION 

Using Figs. [2] and [8] as "reference plots", we may now 
examine the effect of a weak, symmetry breaking ran- 
dom potential V(6). This potential is chosen to consist 
of ten steps of equal widths with a height that is a (uni- 
formly distributed) random number and varies between 
and 0.01. Figure [T4l shows the specific random poten- 
tial chosen. In this case we start with perfectly symmet- 
ric/antisymmetric configurations, a = ±1. 

The time evolution of the symmetric configuration 
shown in Figs. 2] and [5] is almost identical with that of 
Figs. [2] and [3j i.e., the case considered in the previous 
section with V — 0. The two blobs merge rather rapidly. 
The antisymmetric case shown in Figs.[TU] and [11] is of 
greater interest. Here, after a relatively short time, the 
system passes through a "quasi-equilibrium" configura- 
tion, seen as the plateau in the plot of energy versus 
time in Fig. 1111 During this time interval, there are two 
localized blobs. However, parity is no longer a conserved 



quantity in this case. There is no symmetry in the system 
to preserve the node that was built into the initial con- 
ditions. As a result, the two blobs eventually merge into 
one in contrast to the results of Fig. [51 In other words, 
the apparent repulsion of the two solitary waves is not 
present at sufficiently large t. 



V. EFFECT OF SLIGHT ASYMMETRIES IN 
THE INITIAL CONFIGURATION 

In another set of runs, we set the random potential to 
zero and select slightly asymmetric initial configuration, 
with a = ±1.01. Our initial condition is thus not a par- 
ity eigenstate and our calculations show that again the 
two initially distinct blobs merge after a characteristic 
timescale. The qualitative features of this calculation are 
the same as in the case of a random potential described 
in the previous section. 

In the case of an almost symmetric initial configura- 
tion, a — 1.01 shown in Fig. [SI the two separate blobs 
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FIG. 12: Snapshots of \^f(9, t)\ (higher curve in each panel) 
and of 4>{9,t) (lower curve in each panel) of the order param- 
eter * = \^(8,t)\e i<t ' { - B ' t \ for the almost antisymmetric initial 
condition, a = -1.01, for t/to = 0, 10, 100, 300, and 400. The 
axes are the same as in Fig. 1. In all the above graphs there 
is no external potential, V = 0. 



ative) parity eigenstate, and that dissipation must have 
been weak. 

It is instructive to estimate the timescale, to, for our 
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FIG. 14: The random potential in units of E = h 2 /(2MR 2 ) 
as function of the angle 6 that is used in the specific calcula- 
tions shown in Fies.lH and 1101 
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FIG. 13: The energy of the system in units of Eq — 
h 2 /(2MR 2 ), versus time for < t/to < 400 corresponding 
to Fig.Hl 

merge rapidly, very much as in Figs. [2] and [4] On the 
other hand, the almost antisymmetric case, a = —1.01 
shown in Figs.fT^l and [T51 shows a plateau in the energy 
and a period of "quasi-equilibrium" during which the two 
blobs have relatively well-determined shape and location. 
Eventually, however, the two blobs merge again into a 
single solitary wave, much as in Fig. [10] but unlike Fig. [8] 

VI. DISCUSSION AND CONCLUSIONS 

According to the results of our study, the observed re- 
pulsion between bright solitary waves in the experiment 
of Ref. [f| implies that the conservation laws were not 
substantially violated during the time interval investi- 
gated. More precisely, it suggests that deviations from 
axial symmetry in the trapping potential must have been 
small, the initial configuration was very close to a (neg- 



study. If one considers a value of R equal to the lon- 
gitudinal size of an elongated trap, R ~ 0.1 mm, then 
to ~ 10 sec, which is a rather long timescale for these ex- 
periments. Therefore, it seems likely that the character- 
istic timescale over which the experiment of Ref. [5] was 
performed was significantly smaller than the timescale re- 
quired to see the separate blobs merge. Higher tempera- 
tures would enhance the dissipation in the gas and would 
decrease the characteristic time that is required for the 
blobs to merge. To the extent that the deviations from 
axial symmetry in the trapping potential and the anti- 
symmetry in the initial configuration considered here are 
representative of the actual experimental situation, our 
results support the explanation offered in Ref. Q ■ Direct 
experimental determinations of these quantities and of 
the strength of dissipation would thus be welcome. It 
would also be of interest to investigate the long time sta- 
bility (or instability) of the configurations observed in 
Ref. 0. 

The questions examined here may also have important 
consequences on possible technological applications. For 
example, propagation of such solitary waves in waveg- 
uides may serve as signals that transfer energy or infor- 
mation. Therefore, understanding and possibly control- 
ling the way that such waves interact with each other may 
be important. Recent experimental progress in building 
quasi-one-dimensional and toroidal traps should make 
such experiments easier to perform and worth investi- 
gating. 
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